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ABSTRACT 

We present numerical studies of 3-dimensional magnetohydrodynamic (MHD) turbulence in a 
strongly magnetized medium in the extremely relativistic limit, in which the inertia of the charge 
carriers can be neglected. We have focused on strong Alfvenic turbulence in the limit. We have 
found the following results. First, the energy spectrum is consistent with a Kolmogorov spectrum: 

E(k) ~ fc~ 5 / 3 . Second, turbulence shows a Goldreich-Sridhar type anisotropy: k\\ oc k 2 ^ 3 , where 
fen and k± are wavenumbers along and perpendicular to the local mean magnetic field directions, 
respectively. These scalings are in agreement with earlier theoretical predictions by Thompson & 
Blaes. 

Subject headings: MHD — turbulence — relativity 



1. INTRODUCTION 

In a strongly magnetized extremely relativistic 
medium, as in the magnetospheres of pulsars and black 
holes, energy density of the conducting matter can be 
neglected. When this is the case, the Lorentz force van- 
ishes and, thus, we can use force-free approximation for 
the medium. Turbulence in this force-free magnctohy- 
drodynamics (MHD) limit is the topic of this paper. The 
magnetospheres of pulsars and black holes (see Goldre- 
ich & Julian 1969; Blandford & Znajek 1977; Duncan 
& Thompson 1992) or gamma-ray bursts (see Thomp- 
son 1994; Lyutikov & Blandford 2003) are examples of 
relevant astrophysical systems. 

The force-free MHD admits two normal modes, Alfven 
and fast modes (see, for example, Thompson & Blaes 
1998; Komissarov 2002). In the presence of a strong 
mean magnetic field Bo = Bqz, nonlinear interactions of 
these modes result in turbulence cascade. In this paper, 
we focus on Alfvenic turbulence. 

In the non-relativistic limit, Goldreich & Sridhar 
(1995) model provides a good theoretical description of 
strong Alfvenic turbulence. The model predicts a Kol- 
mogorov energy spectrum (E(k) oc fc~ 5 / 3 ) and a scale- 
dependent anisotropy (fc|| oc k 2 / 3 ). Here fc|| and k± 
are wavenumbers along and perpendicular to the local 
mean magnetic field directions, respectively. The scale- 
dependent anisotropy states that anisotropy is more pro- 
nounced at smaller scales. The model was first nu- 
merically confirmed by Cho & Vishniac (2000). Subse- 
quent numerical studies and further discussions are given 
in Maron & Goldreich (2001), Cho, Lazarian, Vishniac 
(2002), Cho, Lazarian, Vishniac (2003). 

Thompson & Blaes (1998) theoretically studied physics 
of force-free MHD turbulence. They introduced two new 
formulations for the force-free MHD and discussed var- 
ious aspects of turbulence processes, including those of 
Alfvenic turbulence. They found that force-free Alfvenic 
MHD turbulence exhibits scaling relations very similar 
to the non-relativistic ones. That is to say, regarding 
to strong force-free Alfvenic MHD turbulence, they ar- 
gued that it has a Kolmogorov spectrum and a Goldreich- 
Sridhar type anisotropy. 



In this paper, we numerically study strong force-free 
Alfvenic turbulence in flat space. In §2, we describe our 
numerical methods. In §3, we present results. We give 
discussion in §4 and conclusion in §5. 

2. METHOD 

We solve the following system of equations: 



where 





Q — (Si, #2, S3, B 2 , B 3 ), 
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S = E x B, 
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E = — l —S x B, 
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where E is the electric field and S the Poynting flux vec- 
tor (see Komissarov 2002). One can derive this system 
of equations from 

dfj* F^ u = (Maxwell's eq.), (7) 

Q^U = _JU ( Maxwell >g eq /) ; ( 8 ) 

dpT^jt = (energy-momentum eq.), (9) 

FvnU 1 * — (perfect conductivity), (10) 

where *F tlL/ is the dual tensor of the electromagnetic 
field, the fluid four velocity, and TPZ the stress-energy 
tensor of the electromagnetic field: 

T^=F^F av -^(F a0 F^)g^, (11) 

where g^ is the metric tensor and F a P is the electro- 
magnetic field tensor. We ignore the stress-energy tensor 
of matter. We use fiat geometry and units with c — 1. 
Greek indices run from 1 to 4. One can obtain the force- 
free condition from Maxwell's equations and the energy- 
momentum equation: d^T^ = —F VI1 J^ = 0. From 

eauation H10(l . one can derive 

E-B = 0, (12) 

B 2 - E 2 > 0. (13) 
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We solve equations (l)-(6) using a MUSCL-type 
scheme with HLL fluxes (Harten et al. 1983; in fact 
in force-free MHD these fluxes reduce to Lax-Friedrichs 
fluxes) and monotonized central limiter (see Kurganov 
et al. 2001). The overall scheme is second-order accu- 
rate. After updating the system of equations along x 1 
direction, we repeat similar procedures for x 2 and x 3 
directions with appropriate rotation of indexes. Gam- 
mie, McKinney, & Toth (2003) used a similar scheme for 
general relativistic MHD and Del Zanna, Bucciantini, & 
Londrillo (2003) used a similar scheme to construct a 
higher-order scheme for special relativistic MHD. 

While the magnetic field consists of the uniform back- 
ground field and a fluctuating field, B = Bo + b, electric 
field has only fluctuating one. The strength of the uni- 
form background field, Bq, is set to 1. We use a grid of 
512 3 . At t = 0, only Alfven modes are present in the 
range 



4 < k± < 6 and 
1 < k\\ < 2 



(14) 
(15) 



in wavevector (k) space. The MHD condition E • B = 
is satisfied at t = 0. The energy density of the random 
magnetic and electric fields at t = is ~ 0.1. Therefore, 
we have 

- bkl - 1 Met\ 
X = s-, 1 (16) 



Boh 



at t = 0. 



3. RESULTS FOR STRONG TURBULENCE 

Figure 1 shows energy spectra of magnetic field. At 
t = (not shown) only large scale (i.e. small k) Fourier 
modes are excited. At later times, energy cascades down 
to small scale (i.e. large k) modes. After t ~ 3, the 
energy spectrum decreases without changing its slope. 
The spectrum at this stage is very close to a Kolmogorov 
spectrum: 

E(k) oc k~ 5/3 . (17) 

In Figure 2, we plot contour diagram of the second- 
order structure function for magnetic field in a local 
frame, which is aligned with the local mean magnetic 
field B L : 



SF 2 (r 



II > r ^ 



=< |B(x + r) -B(x)| 2 > avg 



(18) 



where r = riirii + rj_fj_ and f|| and fj_ are unit vectors 
parallel and perpendicular to the local mean field Bl, 
respectively. See Cho et al. (2002) and Cho & Vishniac 
(2000) for the detailed discussion of the local frame. The 
contour plot clearly shows existence of scale-dependent 
anisotropy: smaller eddies are more elongated. By ana- 
lyzing the relation between the semi-major axis (~ l\\ ~ 
1/A;||) and the semi- minor axis (~ l± ~ l/fc±) of the 
contours, we can obtain the relation between fc|| and k±. 
The result in Figure 3 is consistent with the Goldreich- 
Sridhar type anisotropy: 

fc|| cx k 2 / 3 . (19) 

The electric field shows similar scalings as the magnetic 
field. Our numerical results for spectrum and anisotropy 
are consistent with theoretical predictions by Thompson 
& Blaes (1998). 



4. DISCUSSION 
4.1. Phenomenology 

Thompson & Blaes (1998) derived scaling relations for 
force-free MHD turbulence, which are confirmed by our 
simulations. Here we re-derive the scaling relations using 
a simple phenomenology. 

The magnetic and the electric fields follow the equa- 
tions: 

— = -V x E, (20) 
<9E 

w VxB •'• 

where the current density J is 

_ (E x B)V - E+(B-VxB-E-Vx E)B 

J ~ & 
(Lyutikov 2003). The current density J in equation 121fl 
is essential for nonlinear interactions. In equation 1)22(1 . 
when Bq 3> b and k±^ 3> fcii , the first term on the right 
side dominates for Alfven turbulence. 

Suppose that we have a Alfven wave packet whose par- 
allel s 



(21) 



(22) 



is 2 1| ~ ku 1 and perpendicular size l± 

(~ / ~ when anisotropy is present). This wave 
packet travels along the magnetic field line at the speed of 
c (= 1). When this wave packet collides with opposite- 
traveling wave packets (of similar size), the change of 
energy (AS) per collision is 

AS ~ (dS/dt)At ~ (kbf/B Q )/k h (23) 

where we use equation ((21f) to estimate d£/dt. We as- 
sume Bq ^> b and k± ^> kn, and, hence, J ~ kb 2 / Bq (see 
equation iJUJl). Note that b\ ~ J5j. Therefore, 

AS/S ~ (kbi)/(k\\Bo) - t w /t eddy , (24) 
where t w (= l\\/c— 1/fcy) is the wave period and t e dd y (= 
l±/vi ~ l/(k±b/B )) the eddy turnover time. Here we 
used vi ~ cb/B = b/B . We consider two cases: x = 
(kb^/ihBo) - 1 and X < 1. 

If x (= kbi/k\\Bo) ~ 1, the resulting turbulence is 
strong (see Goldreich & Sridhar 1995 for more discus- 
sion about non-relativistic Alfvenic turbulence and the 
role of x)- O ur simulation satisfies this condition (Fig- 
ure 4). Since x ~ tw/teddy (equation |24j ) . the condi- 
tion x ~ 1 implies that there is a balance between the 
hydrodynamic and the wave times-scales, which means 
that we can use either time-scale for the energy cascade 
time-scale. Therefore, the constancy of spectral energy 
cascade rate in this case becomes 
b 2 

— l — ~ bfvik cx bfk = constant, 

teas 

which results in 

bi oc l 1/3 or E{k) oc fc~ 5 / 3 , 

where we use kE(k) 
obtain 



bf. From kbi/(k\\B ) 



(25) 

(26) 
1, we can 



,2/3 



-II - » ± , (27) 
where we use k ~ fcj^. 

If X < 1) then the turbulence is weak: many colli- 
sions are required to make A£/£ ~ 1. Thompson & 
Blaes (1998) argued that force-free MHD turbulence in 
this regime has virtually constant and E(k) oc fc -2 . 
See Galtier et al. (2000) for non-relativistic weak MHD 
turbulence. 
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Fig. 1. — Spectra at t=3.15 (dotted) and t=5.5 (solid). The spectra are compatible with a Kolmogorov spectrum: E(k) oc k 5 / 3 . 

Fig. 2. — Contour plot of the second structure function for the magnetic field at t = 4.7. Contours, representing eddy shapes, clearly 
show scale-dependent anisotropy: smaller contours are more elongated. 




SEMIMIN. AXIS («1 /k | ) 



Fig. 3. — Anisotropy. Semi-major and semi-minor axes are obtained from Figure 2. We observe a Goldreich-Sridhar type anisotropy: 
fcii oc k 2 ^ 3 . 

Fig. 4.— x (= kbi/k^Bo ~ ll\\/l±][SF 2 (0,l\\)] 1 / 2 ) is nearly constant in our simulations. We take the ratio l\\/l± from contours in Figure 
2. Note that the x-axis is ix(~ V^J-) m S r id units. 



4.2. Scaling of fast modes 

At the beginning of the simulation, we have only Alfvcn 
modes. As later times nonlinear interactions of Alfven 
wave packets produce fast modes. Figure 5 shows the 
amount of fast modes generated from the Alfvcnic tur- 
bulence. The ratio of fast to Alfven energy is roughly 
0.13 to 0.15, which can be a measure of mode cou- 
pling between Alfven and fast modes. See Thompson 
& Blaes (1998) for detailed discussion of mode cou- 
pling between Alfven and fast modes. It is interest- 
ing that the energy ratio is not very much different 
from the energy ratio of compressible to incompressible 
modes in hydrodynamic (Porter, Woodward, & Pouquet 
1998), super- Alfvenic non-relativistic MHD (Boldyrev, 



Nordlund, Padoan, 2002), or mildly sub- Alfvenic non- 
relativistic MHD (Cho & Lazarian 2002) cases. Un- 
like our current simulation, all these works used either 
isotropic initial condition or isotropic driving. 

Cho & Lazarian (2002) showed that fast modes gen- 
erated from solenoidal components in non-relativistic 
MHD turbulence are isotopic when driving or initial 
condition is isotropic. In case of anisotropic driving 
or an anisotropic initial condition, fast modes can be 
anisotropic. But we expect that anisotropy of fast modes 
in that case is scale- independent. 

Can we derive a similar conclusion for force-free MHD? 
In Figure 6, we plot contour diagram of the second-order 
structure function for fast modes. As in the case of 
Alfven modes, we can assume that contours represent 
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Fig. 5. — Generation of fast modes. Energy in fast modes grows rapidly initially. After t > 0.5, the energy ratio of fast to Alfven modes 
is about ~ 0.13. 

Fig. 6. — Contour diagram of the second-order structure function for fast modes at t = 4.7. Contours show anisotropy that is virtually 
scale-independent . 



eddy shapes. Contours do show anisotropy. But, 
anisotropy does not show strong scale-dependence. 



the 



4.3. scaling of current 

In equation l|21|l. we are not interested in the parallel 
components of the right hand side of the equation. Per- 
pendicular components are important for the evolution 
of the electric field. Therefore let us consider only the 
perpendicular components of the equation. We also as- 
sume that we have only Alfven modes. Then the three 
terms on the right hand side of equation ill' 2 1) can be 
approximated by 



(E x B)V • E/B 2 ~ kbf/B , 
(B • V x B)B/B 2 ~ k\\b 2 /B , 
(E - V x E)B/B 2 ~ kbf/B 2 , 



(28) 
(29) 
(30) 



where we assumed b\ ~ E\. Therefore, the dominant 
term is the first term. From this we can easily show that 
the spectrum of the current density follows 



Ji oc r 1 / 3 , 



or 



Ej(k) oc k- 1/3 . 



(31) 



On the other hand, the charge density follows 

p e oc V • E oc fc 2/3 , (32) 

or 

E Pte (k) oc fc 1/3 . (33) 

5. CONCLUSION 

Using numerical simulations, we have studied 3 dimen- 
sional force-free MHD turbulence. We have found that 
energy spectrum is compatible with a Kolmogorov spec- 
trum: E(k) oc fc~ 5 / 3 . We have calculated anisotropy and 

found a Goldreich-Sridhar type anisotropy: fcii oc k 2 / 3 . 
These findings are consistent with earlier theoretical 
studies in Thompson & Blaes (1998). 



I thank Chris Thompson for useful discussions and ad- 
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sions. This work utilized CITA supercomputing facili- 
ties. 
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